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Fundamental rules and definitions of Fractional Differintegrals are outlined. Factorizing 
1-D and 2-D Helmholtz equations four semi-differential eigenfunctions are determined. The 
functions exhibit incident and reflected plane waves as well as diffracted incident and reflected 
waves on the half-plane edge. They allow to construct the Sommerfeld half-plane diffraction 
solutions. Parabolic- Wave Equation (PWE, Leontovich-Fock) for paraxial propagation is fac- 
torized and differential fractional solutions of Fresnel-integral type are determined. We arrived 
at two solutions, which are the mothers of known and new solutions. 



1 Introduction 

The mathematical theory of the fractional calculus and the theory of ordinary frac- 
tional differential equations is well developed and there is a vast literature on the subject 
[1] , [2] , [3] , [4] , [5] and [6] . The theory of partial fractional differential equations is a 
recently investigated problem and the theory mainly concerns fractional diffusion-wave 
equations [7], [8], [9], [10] and [11]. 

The main objectives of this paper is a factorization of the Helmholtz equation to ob- 
tain four semidifferential eigenfunctions allowing to construct the well known half-plain 
diffraction problem. Factorizing the Leontovich-Fock equation, we determine semidiffer- 
ential Green functions, which allow to find paraxial solutions for a given beam boundary 
conditions. 

The article is organized as follows. In Sec. 2, we quote the required rules for fractional 
differintegrals and four fundamental definitions of fractionals. Section 3 is devoted 
to a factorization of ordinary differential equations and determinations of fractional 
eigenfunctions. Section 4 and 5 constitute the main body of our paper and contain 
derivations of fractional eigenmodes for the Helmholtz equation and fractional solutions 
of Leontovich-Fock equation. The final section is devoted to comments and conclusions. 



Z Main Rules and Definitions of Fractional Differintegrals 

Integration and differentiation to an arbitrary order named fractional calculus has a 
long history, see[l], [3]. In the background of the fractional operations, we see general- 
ization of integer order calculus to a noninteger order, class of differintegrable functions 
and applications of the calculus. We recall required rules of fractional diffcrintegral 
calculus. For simplicity, at present let us assume that the considered functions are real 
and differintegrals are of real order and are related ti the interval of real axis from c to 

CO. 

1. Analyticity: If functions / : R — ► R is analytic for x € (c, oo) then the 
fractional of g-order prescribed to the interval (c, oo) ; c D q z (f(x)) is an analytic function 
of x and q. 
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2. Consistency: Operation c D^.(f(x)) must be consistent with the integer order 
differentiations if q = n and with integer order integrations if q = —n. The operation 
must vanish together with n-1 derivatives at point x = c. 

3. Zero Operation: c D°(f(x)) = f(x) 

4. Linearity: c D%(af(x) + bg(x)) = a c D%.(f(x)) + b c D%{g(x)) where a and b are 
arbitrary constants. 

5. Additivity and Translation: c D%( c Dl(f(x))) = c 

The definition of fractional diffcrintegrals due to Griinwald and Letnikov is the most 
fundamental in that it involves the fewest restrictions on the functions to which it is 
applied and avoids the explicit use of notations of the ordinary derivative and integral. 
The definition unifies two notions, which are usually presented separately in classical 
analysis; the derivative and integral (differintegral) . It is valid for a q-th order derivative 
or (-g)-th folded integrals whether or not q is a real integer. Thus the differintegral of 
f:R — >R is: 

h-q N_1 re ' — i 

((2.D) Mm) = „Hm ^ g J^f )f{ x - jh), 

where h = (x — c)/h, and x > c. The definition stems from the difference quotient 
defining ra-th order derivative, which contains n terms in nominator and n-th power of 
h in denominator. The number of terms tends to infinity for noninteger q in the nomina- 
tor and g-th power of h in the denominator. Just like in the case of binomial formula for 
the positive integer power (a + b) n and for negative power as well as noninteger power 
(a + b) q . The convergence is a critical point but the formula (2.1) is very convenient for 
computations. The definition is equivalent to the Riemann-Liouville fractional integral: 

(JVM*) = c D-"(f(x)) = ^J\x-ty- 1 f(t)dt, ((2.2)) 
< q < 1 
and fractional derivative: 

Mm) - r^ljy*-^ 9 ^- « 2 - 2 » 

< q<l. 

The definitions stem from consideration of Cauchy's integral formula and are very 
convenient to implement but require observation of convergence of the integrals. Nev- 
ertheless, because of their convenient formulations in terms of a single integrations they 
enjoy great popularity as working definitions. They play an important role in the devel- 
opment of the theory of diffcrintegrals and for their applications in pure mathematics- 
solutions of noninteger order differential equations, definitions of new function classes, 
summation of series, etc. In this paper, we do apply the these definitions in spite of 
some ambiguity in respect to consistency and additivity rules. However, the demands of 
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modern applications require definitions of fractional derivatives allowing the utilization 
of physically interpretable initial conditions, which contain classical initial conditions for 
functions and their integer derivatives at the initial point. Unfortunately, the Riemann- 
Liouville approach leads to initial conditions containing fractional derivatives at the 
lower limit of integrals. A certain solution to this conflict was proposed by M. Caputo 
[4], whose definition is as follows: 



t)dt, ((2.3)) 



< q<l. 



where C denotes the Caputo definition. 

For the purpose of applications to fractional differential equations, we introduce the 
Miller-Ross sequential fractional derivatives. The main idea is based on the relation: 



D na f(x)= p a D a ...D a J (x) 

n 

or D a f(x) = D ai D a2 ...D a "f{x) 
where a = a.\ + 012 + ••• + a„ 

and the simplest fractional equation of order 0(N, q) takes the form: 



N 

((2.4)) = /(*>») , 

where v = -and the adequate initial conditions for fractional derivatives, see [3] and 
[4]- 

We shall call (2.4) the fractional linear differential equation with constant coeffi- 
cients of the order (N,q), where q is the least common multiple of the denominators 
of the nonzero aj = jv. The solution to the equation can be found by use of Laplace 
transformations. We know that N-th order linear differential equation has iV linearly 
independent solutions. In [3], it is shown how to construct linearly independent Ab- 
solutions of homogeneous fractional differential equations. 

3 Factorization and Eigenfunctions of an ODE 

An eigenfunction y(£) of the linear operator L[y(^)] is such a function that the repeated 
operations preserve the function, e.g. L[y(£)] = Cy with the exactness to a multi- 
plicative constant C . In the case of fractional operation, the definition is extended 
to preservation of the function but an additive constant or a term of power of £, e.g. 
(7r£) -1 / 2 , is subtracted at each step, see [1] and [4]. 
Consider the following ODE: 
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((3-1)) DZV(0=V(0> 

which possesses two eigenfunction y — e ± ? and two eigenvalues ±1 . Factorizing the 
last equation according to: 



(a 4 -l) = (a 2 + l)(a 2 -l), 
a 2 -I = (a+l)(a-l) 

and 

a 2 + 1 = (a + i)(a - i) 
where a = D 1 / 2 , the solution to the semi-differential equation, takes the form: 

tf(0 = Ito(0 + £» 1/2 l/o(0. ((3-2)) 
2/o (0 = e« 

and y(£) satisfies the semi-differential equation; D 1 / 2 y — y = as well as the equation; 
-D 2 J/ ~ 2/ = according to the above mentioned rules, that is consistency and additivity. 
On the other hand, for y = , we have the fractional eigenmode: 

which satisfies the equations; D x l 2 y + iy = and D 2 y + y = 0. According to the 
Riemann-Liouville definition the semi-derivative on the interval x G (0, oo) is given by 
the formula; 

otfJ/V = -^Bi,i/2(0 = -^ + e^Erf{^l), 
where E 11 / 2 (£,) is a Mittag-Lcfflcr function, sec [4]. The eigenfunction of the equation 

1/2 

(oD^ y — y = 0) takes the form; 

((3-3)) y(0 = e« + ^i,i/ 2 (0 = ^ + e*Erfc{-y/i), 

where the complementary error function is: 

oo 

Erfc(-y/Z) = ^= J eM~t 2 )dt. 

This is a well known solution, see [4], and it is an eigenfunction of the operator (3,1) in 
the sense of the extended definition of fractional eigenfunctions. The problem is that 
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by substitution of (3.3), we have: 

-!>»<« = "v^. 

In the next Section, we discuss the deficiency and remove power terms in the case of 
the 2-D Hclmholtz equation. 

4 Factorization and Eigenfunctions of Helmholtz Equa- 
tions 

Let us consider the following 1-D Hclmholtz equation: 
((4.1)) D 2 x y(x) + k 2 y(x) = 0. 

OO 

Taking x = ikt; we are constructing 4-eigenfunctions by use of Fresnel integrals J exp(it 2 )dt 

X 

instead of the complementary error functions. We note, that in principle, there are the 
following ordinary modes; e lkx , e~ lkx and fractional modes: 

oo- oo 

((4.2)) e lkx J e~ u2 dt, e' lkx J e u * dt. 

yjkx y/kx 

The first and the second pairs of modes are complex conjugate and since there are four 
resulting modes: 

oo oo 

sin(M,cos(M, j cos{kx ^ )dt , J Mkx - t i )dt . 

Vkx \f~kx 

Substituting the two last terms into (4,1), the following terms; —k 2 /2^/kx , ^fkx/Ax 2 
appear on the right hand side of the equation, respectively. The power terms are the 
consequences of the accepted definition of fractal derivatives and in the considered 
case the consistency and additivity rules are not completely preserved. We withhold a 
discussion of any interpretation of the derived fractional modes. But for two-dimensional 
problems we demonstrate a method to remove the ambiguity as the problems have a 
pronounce physical meaning. In the case of the 2-D Helmholtz equation: 

((4.3)) A x ^(x,y) + kH(x,y)=Q, 



it is found, that 
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oo 



((4.4)) &(x,y) = e - iax - iby J e tt2 dt, 

u(x,y) 

with k = \Ja 2 + b 2 , satisfies (4.3) if the function u(x, y) obeys the following characteristic 
equations: 



((4.5)) u[(d x u) 2 + (dyu) 2 ] — ad x u + bd y u, and d XfX u + dy iV u = 0. 

Introducing an analytical function; f(z) = u(x, y) + iv(x, y) with possible singularity at 
z = x + iy = 0, we note, that (4.5) is the real part of the following equation: 

r'(z)f(z)f(z) = Kr'(z), 

where n = a + ib, \n\ = k = \J a 2 + b 2 . The equation can be reduced to, called by us, 
an eikonal equation of diffraction; 

((4.6)) ± {f{z) f = 2n. 

Neglecting the constants of integration, the solution of the last equation is 

f(z) = ±V2kz 



((4.7)) and u(x,y) = ±y/kr + ax + by = ±V2fcrcos( ^ - — ), 

where r 2 = x 2 + y 2 , (r, 9) are polar coordinates and a is the incident angle of a plane 
wave, see [12]. In the case of reflected waves, we have: 

e - iax+iby F( Ul ). 

By use of the eikonal equation of diffraction, we obtained: 

9 + a, 



U\{x, y) = \Jkr + ax — by = V2kr cos( — - — ). 

It deserve notice, that the obtained complex solution for the 2-D equation; f(z) = 
±\/2kz is analogous to 1-D case, where we have ±Vkx for the lower limit of Fresnel 
integrals. Hence, we have 4-eigenmodes of the Hclmholtz equations: 

g—iax—iby ^ — iax+iby 



e- iax - iby F(u(x,y)) and e - iax+iby F( Ul {x, y)), 
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where F(u) = J°° e 4 * 2 dt. Rearranging complex and complex conjugate terms, we obtain 
the full set of real modes of 2-D Helmholtz equation: 

sm(ax ± by) 7 cos(ax ± by) , ((4-8)) 

oo oo 

J sin(aa; ± by — t 2 )dt and J cos(ax ± by — t 2 )dt. 

^/kr+ax±by ^/kr+ax±by 

One may speculate, that a fractional Laplacian operator: 

A 1 / 4 = n £ 1/4 
x.y U ^x,y 

has fractional eigenfunctions represented by two last terms of (4.8) and the operator 
is related to half-plane (x,y) and exhibits half-plan edge waves. Fractional Laplacian 
operations may concern 2-D surfaces, like circle, ring, holes in plane, etc., and they are 
concerned with respective edge waves. 

In the case of the 3-D Helmholtz equation and for the wave: 

/•OO 

e -iax-iby-icz / e 2t 2 d ^ 
J u(x,y,z) 

we derived the following characteristic equations: 

u[{d x uf + (d y u) 2 + [d z uf] = ad x u + bd y u + cd z u, ((4.9)) 
Au(x,y,z) = 0. 

There is no known solutions to these equations. The well known solution, in the 
diffraction theory, the 3-D case: 

g—iax — iby — icz . 



f 

Ju{ 



i(x,y) 

where a = fccosacos/3, b = sin a cos j3, c — fcsin/3 and r = \J x 2 + y 2 , leads to 
u(x, y) = \Jkr cos [3 + ax + by, which satisfies the 2-D characteristic equations and it is, 
in principle, reduction of a 3-D problem by variable separation to the 2-D problem. We 
note, that equations (4.9) are also obeyed. 

5 Fractional Solutions for Paraxial Propagation 

Transition from the rigorous wave theory based on a 3-D Helmholtz equation: 



((5.1)) V 2 U(x, yz) + k 2 s(x, y, z)U{x, y, z) = 0, 

where e is a slowly changing dielectric permeability, to the transversal diffusion ap- 
proximation leads to the change of the kind of the differential equations and to a new 
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formulation of the boundary value problem. In the contrast to the elliptic Helmholtz 
equation, the Schrodinger type parabolic wave equations (paraxial-wave equation) de- 
scribes the evolution of the wave amplitude in process of almost unidirectional propaga- 
tion along the optical axis. Physically, it means neglecting the backward reflections from 
the interfaces and an inaccurate description of the waves diffracted into large off-axis 
angles. This approximate approach finds a wide spectrum of applications in radio wave 
propagation, underwater acoustics, linear laser beam propagation and X-ray imaging 
optics, see [13], [14], [15], [16] and [17]. A variety of its modification has been used in 
the diffraction theory, nonlinear optics and plasmas, see [18], [19], [20], [21], [22] and 
[23]. 

Substituting the following form of expected solutions: 

U(x, y, z) = u(x, y, z) exp(ik * r) 



where k = (p, q,j), |k| = k, -f 2 = k 2 - (p 2 + q 2 ), k * r =px + qy + -fz, sin/3 = fc +<? , 
to equation (5.1) and assuming e = 1 + ot(x,y, z), \a\ <C 1 , sin/3 « , d ZtZ u « 0, we 
derive: 



((5.2)) 2ikd z u + A XtV u + k 2 a(x, y, z)u = 0. 

For the parabolic equation, the Cauchy problem with the given initial distribution 
u(x,y,0) = Uo(x,y) (named also one-point boundary value) is correctly posed if some 
radiation condition is added excluding spurious waves coming from infinity. 

Considering the two cases; (2+l)D and (1+1)D, we neglect the inhomogcncity of 
the dielectric permeability a « 0, and write the two following equations: 



((5.3)) 2ikd z u + A x . y u = 0. 

((5.4)) 2ikd z u + d x . x u = 0. 

Next, we factorize the equation(5.4) to obtain: 

Q D 2 X + 2ik Dl = [o D x + Vk{l - i) DlJ 2 ] * [ Q Dl - Vk{\ - i) Q D]/% 

where the fractional derivative, according to Riemann-Liouvillc definition, takes the 
form 



z 

D l J 2 u{x,z) = -^=<9 Z J ' u(x,£) 



and the formula is used to describe nonlocal surface admittance, see[17]. 

Applying the Laplace transformation with respect to the variable z to the following 
equation: 
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((5.5)) [ Dl ± Vk(l - i) Q D\l 2 ]u{x, z) = 0, 

we obtain: 

((5.6)) [ Dl ± - i) Dl /2 }u(x, s) = F(x). 

and neglecting the initial condition F(x) at z = 0, we find: 

M(x ' z) = v^P exp( ^ } - 

We now factorize (5.6) again and determine the following solution: 



((5.7)) u{x,y) = \\- f cxp(it 2 )dt 

k V 7r 7 



It is easy to check that (5.7) satisfies (5.4) and its first as well as higher derivatives with 
respect to x also satisfy (5.4). Therefore, we call it the mother of solutions. By variable 
separation, we can write the mother of solutions for (2+l)D equation (5.3) in the form: 

oc oo 

((5.8)) u(x,y,z) = -( J exp(it 2 )dt)( J exp(it 2 )dt). 

which is related to fractional Laplacian with respect to a certain region of plane (x, y). 
It will cause no confusion if we use the same notation u to designate different solutions 
of (5.3) and (5.4). The derivatives 9™ 9™ of (5.8), where m and n are natural numbers, 
satisfy the PWE and may be related to higher order Gaussian-Hermite optical beams. 

Let us calculate the derivative d x d y of (5.8) to derive classical paraxial Green func- 
tion: 

It is well known that, in the problem of diffraction by plane screens (e.g. a thin 
zone plate) , the parabolic approximation is equivalent to the simplified Frcsncl-Kirchoff 
diffraction theory. In fact, any solution of (5.3) can be expressed for z > in terms of 
its one-point boundary value (initial value) over an aperture. Assuming the one-point 
boundary condition over the plane (x,y): 

x 2 + y 2 
uo(x, y) = exp( g — )' 



where wq is the beam radius at z = and by use of the convolution integral: 



10 



A. J. Turski, B. Atamaniuk and E. Turska 



oo oo 



((5.9)) U(x,y,z)= j { J u (x,y)G(x-x,y-y,z)dx)dy, 



we derive: 



« 5 - 10 » "'••'•* ) -55<JT ? j« , '' < -T5^>' 

and by simple algebraic manipulation, we can write the classical form of the laser beam, 
see [24] and [25]: 

// -.-.w „, .Wo , r 2 _-, ikr 2 . 

((5.11)) U(x,y,z) — — — exp( — expyi tan +- 



w(z) v w{z) 2 ' ly 2R(z)' 



where r 2 = x 2 + y 2 , w(z) = w oy /l + [z/zr) 2 is the beam radius, zr — 7ru>o/A = 
kw^/2 is the Rayleigh length and R(z) = z + z 2 R j ' z is the Raylcigh curvature. Higher 
order modes can be derived by calculation of the following derivatives: 

(-ir(-ird^u(x, y ,z). 

By inspection or simple reasoning one can see that they satisfy (5.3). In view of the 
convolution properties of (5.9), we can differentiate functions Uo(x,y) or G(x,y,z) to 
obtain the same result 

By a symmetry of (5.3), we find that the equation can be reduced to (1+1)D equation 
(5.4) substituting £ = x + y. In virtue of (5.7), we derive the next mother solution: 



((5.12)) u(x,y,z) = ^ J exp{it 2 )dt, 

v(x,y,z) 

where v(x,y,z) = ^~2^-\j\- The solution (5.8) is related to fractional Laplacian with 
the strip symmetry for \x + y\ < const. Higher order modes derived by the following 
differentiation d x and d x _ x — d XiV = d yiV are as follows: 



k k ,ik(x + y) 2 . , k 2 k , , Ak(x + y) 2 . 

and solutions to one-point boundary value problems, e.g. exp[(x + y) 2 /w 2 ], can be 
derived by use of Fresnel-Kirchoff integral (5.9). By symmetry consideration, see [23], 
it seems that the exhibited mother solutions for PWE exhaust all possibilities. 

We now give an example of the Fresnel solution to the following nonhomogeneous 
PWE: 
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((5.13)) 



(2ikd z + A Xty )V{x,y.z) 



\Jx 2 + y 2 



where a is a constant. 



The solution of (5.13) for the homogeneous initial condition; Uo(x,y) = 0, takes the 
very interesting form of the Fresnel beam: 



((5.14)) 



U(x,y,z) 



a j e 

v(x,y,z) 



~dt, 



where v(x,y,z) = J \/ x 2 + y 2 . 



The real and imaginary part of (5.14) as well as the field intensity \U(x, y, z)\ are 
depicted in Fig. 1-3. 



Fig.l 3-D diagrams illustrating the real part of the Fresnel beam Ref/ = 

oo 

a J cos(t 2 )dt, where v(x,y,z) = ^k/2zy/x 2 + y 2 , k — 2, z — 2 and 

v(x,y,z) 

a=l. 
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Fig. 2. 3-D diagrams illustrating the imaginary part of the Fresnel beam ImU = 

oo 

a J sin(t 2 )dt, where v(x, y, z) — y / k/2zy / x 2 + y 2 , k = 2, z = 2 and a = 1. 

t)(a;,£/,z) 



Fig. 3. 3-D diagrams illustrating the field intencity \U\ = 
y/(RcU) 2 + (ImU) 2 of the Fresnel beam. 

6 Comments and Conclusions 

In principle, the number of cigcnfunctions of the differential Laplace operator can be 
arbitrary. According to factorization and a choice of D 1 /" = a, the candidate number 
is n for the Laplacian. In the case of 1-D and 2-D, we derive 4 semi-differential modes. 
The 2-D Laplacian A I>9 possesses two additional eigenfunctions of the form; 

oo oo 

J sin(aa; ±by — t 2 )dt and J cos(ax ±by — t 2 )dt, 

u(x,y) u(x,y) 

where the lower limit of Fresnel integral u(x, y) is to satisfy the set of characteristic 
equations, see (4.5), for the half plane x,y G (0,oo). It is important to note, that the 
modes are related to a half-plane, like fractional derivatives are related to an interval. 
We recall, as an example, that the fractional derivative of order q of the exponential 
function e Xx related to a half-axis is a Mittag-LefHer function but for the whole axis it is 
the same function multiplied by A to power q; \ q e Xx . We do withhold, throughout the 
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paper, from defining fractional Laplacian operator although we determine the eigen- 
function of the operator related to half-plane, strip and circle. There is an expectation, 
that the right definition must be based on an integral of the Fresnel-Kirchoff type and 
relation to the n-dimensional region, where n is a number of independent variables. One 
may speculate, on the ground of Riemann-Liouville definition, that it may be a two- 
folded convolution integral with respect to x and y such that a Laplace transform with 
respect to the variables gives an anticipated results like in the case of single variable 
functions. Also, one may hope, that the Griinwald-Letnikov definition may be extended 
for the fractional operators. A separate paper will be devoted to existence, uniqueness 
and definitions of fractional Laplacian operators. Here considered operators come out 
from factorizations and lead to results, which satisfy the classical Laplacian. It is not a 
necessary condition for the fractional operators. 

We also mention the 3-D Helmholtz equation but for that problem there is no ex- 
pected solution. We do not know the 3-D region, like in the case of half-plane diffraction, 
for which there is an exact solution (eigenfunctions). The solution derived by variable 
separation (a method reducing the 3-D problem to two dimensions) is not applicable to 
the our requirement. 

Consideration of Leontovich-Fock equation is justified not only by its numerous 
applications but to show that the equation is factorable and the fractional equations lead 
to the mother of solutions. The mother is a function of the Fresnel integral satisfying 
PWE, vanishing at z = and her derivatives 9™ 9™ also satisfy PWE. The notion is not 
trivial as there is expectation that the mother solutions and their derivatives generate 
all possible symmetries of on-axis and off-axis beams in the case of homogeneous PWE. 
The last solution (5.14) of nonhomogeneous equation (5.13) is an illustration of paraxial 
beam of the Fresnel type. 
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